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CALCULATION  Q.E  BQDY-JWAVE  P.ULSE  SHAEES.  WITH  ALLOWANCE  FOB 

FREQUENCY-DEPENDENT  Q 


Paul  G.  Richards 


INTRODUCTION 

The  computation  of  body-wave  synthetics,  in  elastic  Earth  models 
that  vary  purely  with  depth,  is  now  routinely  done  by  several  different 
methods.  In  practice,  the  resulting  synthetics  either  look  the  same,  or 
differ  for  reasons  that  are  quite  well  understood  and  which  indicate  which 
methods  are  fallible  in  certain  situations  (Helmberger  and  Burdick,  1979; 
Aki  and  Richards,  1980;  Kennett,  1983;  Richards,  1985;  Chapman  and 
Orcutt,  1985). 

But,  for  computation  of  pulse  shapes  in  anelastic  models,  there  is  a 
less  satisfactory  experience,  in  part  because  only  a  limited  number  of 
plausible  types  of  attenuation  have  been  explored  to  the  extent  of 
incorporation  into  packages  for  the  computation  of  synthetic  seismograms. 
The  most  common  type  of  attenuation  incorporated  into  synthetics  is  that 
characterized  by  constant  Q  (i.e.  a  Q  that  is  frequency-independent,  but 
which  is  allowed  to  vary  with  depth),  which  is  handled  by  methods 
suggested  by  Carpenter  (1967).  But  seismic  data  indicate  that  in  many 
practical  situations  the  assumption  of  a  constant  Q  is  demonstrably  false. 
Unfortunately,  once  the  false  assumption  is  dropped,  it  is  not  yet  clear 
what  should  replace  it. 

Our  recent  work  in  this  area  has  been  based  on  a  method  of  computing 
synthetics  that  is  effective  in  elastic  media  for  situations  where 
geometrical  ray  theory  is  accurate.  These  include  the  interpretation  of 
body-wave  waveform  data  for  earthquakes  and  explosions  at  distances 
from  about  30°  to  85°.  We  first  describe  this  method,  which  would  hardly 
be  needed  (since  there  are  so  many  other  effective  methods  for  evaluating 
the  predictions  of  ray  theory  in  elastic  media  that  smoothly  vary  with 
depth),  were  it  not  for  the  generalization  it  allows  in  the  case  of  anelastic 
media  -  generalizations  that  other  methods  usually  cannot  supply.  Then 
we  show  an  application  of  the  method  to  an  anelastic  Earth  model  where  Q 
is  assumed  to  follow  a  power  law  frequency  dependence.  The  outcome  is  a 
pulse  shape  with  properties  quite  different  from  what  might  be  expected 
if  the  familiar  constant-Q  pulse  shapes  are  surveyed. 
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Our  conclusion  is,  that  a  fundamental  problem  in  the  interpretation  of 
body-wave  broadband  pulse  shapes  -  namely,  the  appropriate  choice  of 
attenuation  model  to  explain  both  spectral  content  and  pulse  width,  and 
how  these  quantities  are  affected  by  propagation  -  is  unresolved.  But,  the 
conventional  model  that  uses  constant  Q  is  flawed. 


RAY  CALCULATIONS  WITHOUT  RAY  TRACING 

An  elementary  use  of  the  relationship  between  travel  time  ( T ), 
distance  ( X )  and  intercept  time  (  x  )  concerns  fitting  a  ray  between  a  fixed 
source  and  receiver,  and  finding  the  geometrical  spreading,  without 
actually  doing  any  ray  tracing. 

Thus,  in  the  frequency  domain,  the  wavefield  for  a  generalized  ray  in 
a  structure  that  varies  solely  with  depth  can  be  written 

W(X0,  (o)  =  J  f(  p,  co)  exp[io>{  Tip,  co)  -  pXip,  co)  +  pXQ  }  ]  dp  ( 1 ) 

where  Tip,  co)  -  pX(p,  co)  is  the  vertical  integral  of  vertical  slowness  along 
the  ray  path, 

T-pX=  x  =  J  yl[l/v2  -  p2]  dz,  (2) 

ray 

in  which  v(z)  is  the  velocity  (complex,  if  there  is  attenuation,  and  also 
frequency-dependent  if  there  is  dispersion);  and  f(  p,  co)  includes  the 
radiation  pattern  as  well  as  relevant  reflection/transmission  coefficients 
for  the  generalized  ray  of  interest.  There  are  minor  differences  of  detail 
between  plane  stratified  and  spherically  stratified  media.  Richards  (1973) 
obtained  examples  of  (1)  even  for  rays  that  have  turning  points,  rather 
than  reflection  from  a  specific  interface. 

In  order  to  obtain  the  ray-theory  approximation  for  an  array  of 
receivers,  consider  the  following  six  steps: 

(i)  One  fixes  <yand  computes  f  and  x{p)  for  a  set  of  equally  spaced  real 
p-values,  in  a  range  of  the  real  ray-parameter  axis  that  is  close  to  where 
one  expects  complex  saddle  points  to  lie  for  the  range  of  distances  (i.e.  X0  - 
values)  at  which  one  wishes  to  know  W.  Let  the  increment  in  p  here  be  Ap. 

(ii)  At  fixed  X0,  one  searches  through  the  array  of  Real[T  -  pX  +  pX0 )  to 
see  at  which  of  the  sampled  points  in  p  it  has  a  minimum  (or,  a  maximum 
for  certain  kinds  of  ray).  Label  this  discrete  point  as  /7jS.  The  complex 
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saddle  must  presumably  lie  near  this  point  of  the  real  p- axis.  If  there  is  no 
attenuation,  the  saddle  will  lie  on  the  real  p-axis  if  there  is  a  real  ray 
between  source  and  receiver. 

(iii)  Find  the  values  (complex,  if  there  is  attenuation)  of  the  three 
constants  T0,p0,  and  DXDP  that  best  fit  the  sampled  phase  factor  according 
to 


T-pX+pX0  =  T0 -ip-p^xDXDP/2  (3) 

in  the  vicinity  of  pjs.  Note  that  if  one  just  uses  pjs  itself,  and  the  points 

just  before  and  just  after  it,  closed  form  expressions  that  in  practice  are 
often  very  accurate  can  easily  be  given  forp0, 70,  and  DXDP.  Aki  and 

Richards  (1980)  introduce  the  function  J(p)  =  T  -  pX  +  pX0  (see  their  page 
423).  In  terms  of  this  sampled  function, 

7js  -1  =  A  Pjs-1>  =  ro  '  ^js  '  Po  >2  *  DXDP  /2  +  APPo-  DXDP  -  ( Ap )2.  DXDP  /  2 

Jjs  =  Apjs)  =T0  -{pjs-p0j2xDXDP/2 

Jjs+\  =  A Pjs+l)  =  T0  -  ( pjs  -  Po  )2  x  DXDP  /  2  -  4p  .  p0 .  DXDP  -  (Ap)2.  DXDP  1 2. 

Therefore,  the  closed  form  expressions  for  the  desired  constants  are 
DXDP  =  -  [7js+i-2/js  +  Jjs4]/(4p)2  , 

Po  =  -  l/js+1 '  ^js  -ll  /  ( 2-  4p  •  DXDP  ) 

and  finally 


To  ^^js  +  (^js-P0)2^^/2.  (4) 

The  point  of  this  exercise,  of  course,  is  that  p0  is  an  estimate  of  the 
position  of  the  complex  saddle;  T0  is  an  estimate  of  the  complex  travel 

time;  and  DXDP  (signifying  dX/dp  at  the  saddle)  is  the  relevant  constant, 
complex  in  general,  needed  to  evaluate  geometrical  spreading. 

(iv)  Interpolate  from  f(  pjs,  <y)  and  f(  pjs+1,  o> )  to  evaluate  f(  p0,  w). 

(v)  Claim  that  the  saddle  point  approximation  to  W  is 

W (X0,(o)  =  f(  p0,  cd)  exp(icdT0 )  x  V  [2tt  +  (icoDXDP)].  (5) 
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Choose  the  next  value  of  X0,  and  loop  back  to  the  second  step  to  do  various 
distances.  Choose  the  next  value  of  frequency  a>,  and  loop  back  to  the  first 
step.  This  loop  is  necessary  only  if  Q  is  frequency-dependent,  and/or  if 
there  is  allowance  for  body-wave  dispersion.  Also,  for  some  types  of 
anelasticity,  it  may  be  possible  to  abbreviate  this  step  and  estimate  more 
directly  the  global  frequency-dependence  of  the  saddle-point 
approximation. 

(vi)  Finally,  one  can  go  to  the  time  domain: 

W(X0,  t)  =  (1/271)  |  W(X0,  (a)  expt-i ft)  t]  d(0.  ( 6 ) 

Note  that  the  whole  effort  is  accomplished  using  real  values  of 
ray-parameter,  at  the  time-consuming  stage  of  tabulating  f  and  x  at 
discrete  p  values.  Often,  the  dependence  on  frequency  in  (5)  is  sufficiently 
simple  that  (6)  can  be  evaluated  explicitly. 

The  method  allows  easily  for  investigation  of  different 
attenuation-dispersion  pairs.  This  is  handled  just  by  insertion  of  some 
appropriate  rule  for  evaluating  v(z,  o),  and  then  recognizing  that  x  =  x(p,co). 
For  example,  in  a  constant  Q  medium  in  which  the  method  of  Carpenter 
(1 967)  is  used  to  handle  dispersion,  the  real  body-wave  velocity  ve(Z)  of  an 
elastic  model  is  replaced  by 

v(z,o»  =  ve(z){  1  +  [  1/Q(g>)]  [(  1  /7t)ln  (  oy<S)0)-  i/2]  }  (7) 

The  method  is  simple  and  rapid  to  execute,  and  has  been  found  quite 
accurate  for  1-D  (that  is,  purely  depth-dependent)  problems  (Richards, 

1985) .  A  comparison  of  synthetics  produced  by  this  method  and  Cagniard 
first  motion  method  is  given  here  in  Figure  1 . 

A  series  of  papers  by  Borcherdt  (most  recently,  Borcherdt  et  al., 

1986)  has  drawn  attention  to  the  need  for  more  careful  handling  of 
attenuation.  Specifically  he  has  advocated  working  with  plane  waves  that 
may  propagate  in  a  direction  that  differs  from  the  direction  of  most  rapid 
attenuation.  (The  "direction  of  propagation"  is  the  direction  of  most  rapid 
phase  increase.)  He  has  shown  that  such  waves  have  properties  that  much 
conventional  analysis  (i.e.,  that  based  on  plane  waves  propagating  in  the 
direction  of  maximum  attenuation)  cannot  reproduce.  Fortunately,  the 
method  based  on  (1)  through  (5),  which  finds  stationary  values  of  the 
ingrand  at  complex  values  of  horizontal  slowness,  gets  around  Borcherdt's 
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valid  objections  to  conventional  analysis,  yet  does  so  for  attenuating 
media  in  a  way  that  requires  little  change  from  computational  experience 
with  elastic -media.  The  subject  is  further  discussed  in  Richards  (1984). 
Incident  lly,  we  note  here  that  the  whole  procedure  has  recently  been 
generalized  so  that  it  carries  over  to  3-D  problems  in  layered  structures, 
with  planar  interfaces  that  can  have  any  strike  and  dip.  In  effect  one  can 
develop  a  generalized  ray  theory  for  a  medium  built  up  from  a  stack  of 
arbitrary  wedges,  and  the  geometrical  ray  approximation  can  be  computed 
by  a  method  that  is  an  extension  of  (1)  to  (7)  (Richards  -  1987  proposal  to 
the  National  Science  Foundation). 


0.6m/i 


Fig.  1 .  Each  pulse  gives  the  P-wave  group  (P+pP  +  sp)  at  distance  80°,  and 
30°  off  the  strike  direction,  from  a  point  source  15  km  deep.  Rake  is  90°, 
and  results  are  given  for  three  different  dips.  These  are  displacements  in  the 
upward  vertical  direction  (after  a  free  surface  correction),  t*  is  taken  as  1 
second,  and  the  time  dependence  of  the  source  is  given  at  the  bottom  of  the 
figure. 

The  upper  row  uses  a  Cagniard  first  motion  approximation  [Langston  and 
Helmberger,  1975].  The  lower  row  (which  has  the  same  time  scale  but  is 
displayed  with  a  slightly  different  amplitude  scale)  uses  the  complex  saddle 
point  spectral  method  reviewed  in  the  text,  with  model  PEM-C  and  a  Q  struc¬ 
ture  that  gives  t‘  =  1  at  80°. 

Numbers  displayed  are  displacement  in  mu.  The  two  methods,  each  of 
which  is  basically  geometrical  ray  theory  adapted  to  an  attenuating  medium, 
agree  to  within  about  10%. 
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EXAMINATION  OF  BODY-WAVE  PULSE  SHAPES  IN  DIFFERENT  TYPES  OF 

MEDIA 


We  have  used  the  method  based  on  (1)  through  (7)  to  compare  pulse 
shapes  in  three  different  Earth  models.  We  have  used  a  fixed 
source-receiver  distance,  one  for  which  geometrical  ray  theory  is  usually 
deemed  adequate  even  for  broad-band  body  wave  pulses,  so  that  attention 
can  be  focussed  on  the  single  issue  of  what  effects  result  from  different 
choices  about  the  frequency-dependence  of  Q. 

Thus,  in  Figure  2  is  shown  the  P-wave  pulse  shape  from  an  impulsive 
source  in  Earth  model  SL8,  for  which  Q  is  constant  and  quite  low  in  value. 
Note  the  charateristic  features  of  this  pulse:  high  frequencies  have  been 
selectively  removed,  and  the  pulse  shape  is  broadened  with  an  asymmetric 
rise  and  fall.  The  attenuation-dispersion  pair  is  based  on  equation  (7). 


Figure  2.  The  P-wave  pulse  at  48°  in  Earth  model  SL8,  for  an 

impulsive  shallow  source.  Constant  Q,  Carpenter  law. 
Tick  marks  are  at  1  s  intervals. 


In  Figure  3  is  shown  the  same  pulse,  but  for  Earth  model  AFL  (i.e.  from 
Archambeau,  Flinn  and  Lambert,  1966).  Q  is  again  constant,  but  is 
significai.Jy  higher  than  the  value  for  SL8. 


Figure  3.  The  P-wave  pulse  shape  at  48°  in  Earth  model  AFL, 
from  a  shallow  impulsive  source.  Tick  marks,  Is. 

Since  both  AFL  and  SL8  are  data-driven  anelastic  Earth  models,  the 
first  based  on  relatively  high  frequency  signals  (about  2  Hz)  and  the 
second  based  on  relatively  low  frequencies  (about  0.005  Hz),  it  is  clearly 
unsatisfactory  that  each  model  is  presented  in  terms  of  a 
frequency-independent  Q  structure.  That  is,  Q  is  constant  and  has  low 
value  in  SL8;  Q  is  constant  but  with  a  high  value  in  AFL.  Various 
high-frequency  phenomena  displayed  in  seismic  data  are  not  present  in 
SL8,  and  low  frequency  phenomena  in  seismic  data  are  not  present  in  AFL. 

It  is  therefore  appropriate  to  consider  some  kind  of  interpolation  for 
Q  across  the  frequency  range  whose  extremes  are  represented  by  the 
frequency  values  typical  of  the  data  bases  underlying  these  two  Earth 
models.  We  have  done  this  by  assuming  a  Q  that  has  a  power-law 
dependence  on  frequency, 


Q(  o>)  -  Q0  (  oil  oj0  ) 


1-s 


(8) 


for  some  choice  of  constant  s.  Q0  is  a  constant,  the  value  of  Q  ( co)  at 
frequency  (oQ . 

We  find  from  causality  rules  and  the  evaluation  of  Hilbert  transforms 


that  the  appropriate  replacement  for  equation  (7)  is  now 


v(z,  (o)  ~  ve(z)  {  1  +  [  1  /  2  Q  (&))  ]  [((  Q  (co)  /  Q  q)  1)  tan  (stc/2)  -  i]}.  (9) 

Consideration  of  the  difference  in  frequency  ranges  over  which  SL8 
and  AFL  were  generated,  and  the  values  of  the  two  separate  but  constant  Q 
values  in  each  model,  leads  one  to  a  value  of  about  0.7  for  the  constant  s 
that  appears  in  equation  (8).  Thus,  Q  has  about  a  0.3  power-law 
dependence  on  frequency,  to  fit  the  two  kinds  of  data  base  underlying 
these  models. 


Figure  4.  The  P-wave  at  48°  in  an  Earth  model  with  power-law 
Q,  having  approximately  the  Q  of  SL8  at  frequencies 
for  which  the  attenuation  in  SL8  was  determined,  and 
the  Q  of  AFL  at  frequencies  for  which  the  attenuation 
of  AFL  was  determined.  Impulsive  shallow  source. 

Tick  marks,  at  Is  intervals. 

The  outcome  of  computing  a  pulse  shape  based  on  this  Q(  co)  might 
reasonably  be  expected  to  be  a  pulse  somewhat  intermediate  between 
those  of  Figures  2  and  3,  since  we  have  merely  found  a  way  to  interpolate 
for  Q  across  a  broad  band  of  frequencies  that  is  representative  of  SL8  at 
one  end,  and  AFL  at  the  other.  But  the  actual  outcome  is  different.  It  is 
shown  in  Figure  4,  and  one  sees  a  pulse  shape  that  has  an  even  shorter  rise 
time  than  does  that  of  Figure  3  for  AFL.  The  reason  turns  out  not  to  be 
hard  to  find.  It  is,  that  at  frequencies  even  higher  than  those  which 
charaterize  the  data  base  for  AFL,  equation  (8)  gives  a  higher  value  of  Q 
than  the  constant  Q  value  of  AFL  itself.  Since  it  is  the  highest  frequencies 
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that  determine  the  rise  time,  this  latter  feature  will,  in  the 
frequency-dependent  Q  model,  be  reflective  of  the  assumptions  built  into 
(8)  concerning  the  very-high  frequency  dependence  of  Q.  In  this  case,  Q  is 
unbounded  as  frequency  rises,  which  presumably  is  why  the  rise  time  is  so 
short  in  Figure  4. 

We  therefore  reach  two  conclusions.  First,  that  a  constant  Q 
attenuation  model  is  unsatisfactory  .  (This  has  long  been  known:  such  a 
model  does  not  explain  observations  across  the  observed  range  of 
frequencies.)  Second,  that  its  replacement  is  not  obvious.  Artificial 
features,  associated  with  attenuation  laws  applied  outside  the  range  of 
frequencies  for  which  there  is  a  good  data  base,  can  in  the  time  domain  be 
drawn  into  the  appearance  of  a  pulse  shape  computed  within  the  frequency 
band  at  which  data  is  acquired.  Another  example  of  this  general  problem 
was  given  by  Richards  (1985)  in  connection  with  an  attenuation  law 
appropriate  for  a  spectrum  of  relaxation  mechanisms. 

A  good  research  plan  to  resolve  this  issue  is  to  see  if  high-quality 
broad-band  pulse  shape  data  can  be  acquired,  for  example  from  RSTN 
body-wave  data  of  deep  earthquakes  observed  teleseismically.  It  may  then 
be  possible  to  obtain  the  attenuation-dispersion  pair  directly  from  the 
data.  We  note  a  first  attempt  in  this  direction  (Choy  and  Cormier,  1986). 
It  is  difficult  to  achieve  results,  because  of  the  need  to  remove  source 
effects,  and  effects  of  heterogeneity  of  structure  at  depth  (e.g.  the 
downgoing  slab,  which  is  always  present  as  a  cause  of  the  earthquake 
activity  itself).  However,  we  are  confident  that  an  observational  program 
will  indeed  result  in  suggestions  for  the  attenuation-dispersion  pair  that 
is  appropriate  in  computation  of  synthetics,  and  that  will  narrow  the 
range  of  current  trade-offs  between  source  and  propagation  phenomena. 
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SEISMICITY  OF  CENTRAL  ASIA 


D.W.  Simpson 


Soviet  Earthquake  Catalogs 


The  Institute  of  Physics  of  the  Earth,  Moscow  issues  an  annual  catalog  of  earthquakes  in  the 
USSR  (Zemletriaseniia  v  SSSR,  hereafter  ZSSSR).  These  annual  volumes  contain  articles  and 
catalogs  for  different  subregions  of  the  USSR,  based  on  seismic  regionalization.  The  largest  of  these 
catalogs  is  for  Central  Asia,  the  active  region  along  the  USSR’s  southern  border. 

The  Central  Asia  region  includes  the  territory  of  the  republics  of  Tadjikistan,  Kazkakhstan, 
Uzbekistan,  Kirghizia  and  Turkmenia  (Figure  1).  The  main  physiographic  divisions  included  are  the 
Northern  and  Southern  Tien  Shan  and  Pamir  Mountains,  Tadjik  Depression,  Kyzul  Kum  desert  and 
Hindu  Kush  (Figure  2b).  The  area  covered  by  the  ZSSSR  catalog  extends  beyond  the  borders  of  the 
USSR  into  Iran,  Afghanistan  and  China  to  the  south.  Coverage  of  these  border  regions  becomes  less 
complete,  however,  because  of  the  lack  of  stations  in  southern  azimuths. 

Comparison  of  ISC  and  ZSSSR  catalogs 

The  increased  coverage  with  regional  and  local  stations  results  in  a  Soviet  catalog  that  has  a  lower 
detection  threshold  and  higher  resolution  than  that  of  teleseismic  catalogs  such  as  the  International 
Seismological  Center  (ISC).  Figure  3  shows  the  annual  and  cumulative  numbers  of  earthquakes 
reported  by  ISC  and  ZSSSR  for  the  period  1964  -  1979  and  maps  derived  from  the  two  catalogs  are 
shown  in  Figures  4  and  5. 

The  availability  of  a  detailed  catalog  for  Central  Asia,  compiled  entirely  from  local  data  sources 
provides  an  interesting  opportunity  to  examine  the  completeness  and  accuracy  of  the  ISC  catalog. 
Unfortunately,  the  ZSSSR  and  ISC  catalogs  are  not  completely  independent,  since  some  of  the  Soviet 
data  are  provided  to  and  used  by  ISC.  Only  a  subset  of  stations  are  reported  to  ISC  however,  so  that 
the  ZSSSR  catalog  can  be  assumed  to  be  more  accurate  and,  because  of  the  increased  density  of 
stations,  it  is  certainly  more  complete. 

We  have  used  an  algorithm  based  on  nearness  in  space  and  time  to  identify  entries  for  common 
events  in  the  two  catalogs.  Of  the  3720  events  reported  by  ISC  for  1964-1978,  3602  are  found  in 
ZSSSR.  Most  of  the  events  reported  by  ISC  but  not  found  in  ZSSSR  are  either  of  low  magnitude  and 
located  near  the  southern  limits  of  the  Soviet  catalog  or  are  poorly  determined  locations  (i.e.  reported 
by  less  than  5  stations).  The  only  well  located  events  above  magnitude  5.0  reported  only  by  ISC 
within  the  area  well  covered  by  ZSSSR  are  events  identified  as  peaceful  nuclear  explosions. 

For  those  events  that  are  common  to  both  catalogs,  we  have  determined  the  differences  in  location 
parameters  (dx  -  longitude,  dy  -  latitude,  dr  -  epicenter  (=V(dx2+dy2),  dz  -  depth,  dt  -  origin  time  and 
dm  -  magnitude). 

Epicentral  location  -  Map  views  of  the  differences  in  location(  dx,  dy),  as  a  function  of  number 
of  stations  used  by  ISC,  are  shown  in  Figure  6.  The  origin  is  the  ZSSSR  location  for  each  event  with 
points  placed  at  the  end  of  the  vector  connecting  the  ZSSSR  and  ISC  locations.  The  difference  in 
epicentral  location  (dr)  is  shown  in  Figure  7.  For  locations  determined  with  less  than  10  stations, 
differences  of  over  200  km  are  observed.  The  difference  in  location  decreases  as  the  number  of 
stations  used  by  ISC  increases  (usually  a  direct  function  of  the  magnitude).  The  differences  in  location 
are  smaller  in  the  EW  direction  and  there  is  an  absolute  bias  to  the  north,  both  most  likely  resulting 
from  the  paucity  of  ISC  reports  from  stations  to  the  north.  For  events  for  which  more  than  50  stations 


are  available  in  the  ISC  solution,  the  standard  deviation  in  the  difference  between  the  ISC  and  ZSSSR 
location  drops  to  about  ±12  km,  with  the  ISC  locations  biased  about  10  km  to  the  ESE  of  those  given 
by  ZSSSR. 

Depth  and  origin  time  -  Figures  8  and  9  show  the  differences  in  depth  (dz)  and  origin  time  (dt)  as 
a  function  of  the  number  of  stations  used  in  the  ISC  location.  While  dr  (Figure  7)  shows  the  gradual 
decrease  with  number  of  stations,  both  dz  and  dt  show  broad  and  biased  (>0)  distributions  even 
beyond  100  stations.The  bias  towards  late  origin  times  and  greater  depths  apparent  in  Figures  8  and  9 
results  from  the  well-known  trade-off  between  these  two  parameters  in  the  earthquake  location 
process.  Many  of  the  catalogs  entries  report  no  depths  or  standard  depths  of  0  km  (ZSSSR)  or  33  km 
(ISC).  Figures  10  and  1 1  remove  these  "unknown"  depths  and  attempt  to  classify  the  remaining  data 
in  terms  of  events  which  we  can  be  confident  are  actually  deep  or  shallow  and  those  in  which  one  of 
the  catalogs  is  in  error.  In  Figure  10  we  assume  that  those  events  for  which  both  ISC  and  ZSSSR 
report  depths  greater  than  60  km  are  actually  deep.  In  this  case  the  distribution  of  dt  is  symmetrical 
about  zero.  We  assume  that  events  that  are  assigned  shallow  depths  by  both  catalogs  (but  not  0  km  in 
ZSSSR  or  33  in  ISC)  are  actually  shallow.  Figures  1 1  and  12  show  that  for  these  shallow  events  there 
is  a  clear  tendency  for  the  ISC  origin  times  to  be  later  (average  of  3  s),  which  is  compensated  for 
(Figure  12)  by  the  ISC  depths  being  deeper. 

Magnitude  -  The  Soviet  classification  of  earthquake  size  is  energy  class,  k,  where  k  is  the  log  of 
energy  release  in  joules.  Rautian  (1960)  gives  the  following  relationship  ,  based  on  an  analysis  of 
determinations  of  both  magnitude  and  energy  class  for  regional  and  local  earthquakes: 

m  =  (k-4)/1.8 
or  k  =  1.8m  +  4 

Figure  13  shows  the  relationship  between  k  as  reported  in  the  ZSSSR  catalog  and  m^  for  those 

events  with  magnitude  reported  by  ISC.  The  relationship  found: 

k=  1.93m +  3.29 

is  similar  to  that  reported  by  Rautian. 


SEISMICITY  ANRJQPQG&AKHY 

Many  of  the  topographic  features  of  Central  Asia  can  be  discerned  in  the  patterns  of  seismicity 
from  the  Soviet  Catalog  in  Figure  2  and  5.  Among  the  more  prominent  features: 

a)  There  is  intense  activity  along  the  southern  front  of  the  northern  and  southern  Tien 
Shan.  This  zone,  the  Gissar  Kokshal  seismic  zone,  is  the  site  of  the  largest  earthquakes  in  Central  Asia 
(Leith  and  Simpson,  1986a)  and  marks  the  southern  edge  of  what  can  be  clearly  identified  as  part  of 
the  "original"  Asian  plate. 

b)  Realtively  low  levels  of  seismicity  are  found  within  the  high  mountainous  areas.  The 
Gissar  and  Kokshal  Ranges  of  the  Tien  Shan  are  flanked  on  the  south  and  north  by  high  seismicity, 
but  the  higher  elevations  are  much  less  active.  Relatively  low  levels  of  activity  in  regions  of  high 
elevaiton  are  also  found  in  the  Talas  Fergana  Ranges  and  the  highest  section  of  the  Pamir. 

c)  Relatively  low  levels  of  seismicity  are  found  within  the  major  intermountain  basins 
and  depressions.  Most  prominent  of  these  is  the  Issyk  Kul  Basin,  where  high  activity  along  the  ranges 
bordering  the  basin  clearly  oudines  the  relatively  low  seismicity  of  the  basin  itself.  The  Fergana  Valley 
is  similar  but  less  obvious.  The  Tadjk  Depression  and  Tarim  Basin  have  high  activity  on  their  northern 
edges,  near  the  Gissar  Kokshal  zone,  but  are  less  active  in  their  interiors. 

The  seismicity  is  thus  concentrated  along  the  major  topographic  gradients  dividing  the  distinct 
physiographic/tectonic  provinces  of  the  region,  with  relatively  low  activity  in  both  high  plateaus  and 
low  basins  and  higher  activity  along  the  flanks  of  the  mountainous  areas. 


Digital  Terrain  Data 


To  study  the  relationship  between  topography  and  seismicity  in  more  detail,  we  have  obtained 
digital  terrain  data  for  Central  Asia  from  the  Defense  Mapping  Agency  (DMA)  topographic  database. 
Topographic  elevations  are  provided  at  3"  intervals  (1200  points  per  degree  or  approximately  100  m 
spacing). 

We  have  carried  out  preliminary  analysis  of  segments  of  these  data  making  extensive  use  of  the 
International  Imaging  Systems  (HS)  image  processor  at  Lamont  Doherty.  The  original  data  are 
provided  in  blocks  of  1°  x  1°  (1200  x  1200  points).  While  that  resolution  has  been  useful  in  some 
detailed  studies  (e.g.  the  Kulyab  salt  dome  study  described  below  or  as  part  of  our  work  on  induced 
seismicity  at  Nurek  reservoir),  the  volume  of  data  is  too  large  for  synoptic  studies  of  large  areas  and  it 
is  necessary  to  severely  decimate  the  data.  The  image  processing  system  provides  a  convenient 
environment  for  reading  the  data,  creating  and  examining  full  resolution  images,  decimating  into 
smaller  panels  and  creating  mosaics  of  larger  map  areas.  The  image  processor  can  also  be  used  to 
apply  different  enhancement  techniques  to  the  topographic  images,  combine  the  topography  with  other 
(eg.  epicenter)  data  and  show  the  results  in  glowing  color.  In  this  report  we  limit  ourselves  to  low 
resolution  greyscale  versions  of  these  images. 

Tadjik  Depression  and  Southern  Tien  Shan 

An  example  of  an  image  created  from  the  terrain  data  is  shown  in  Figure  14.  This  area  includes  the 
Southern  Tien  Shan  (Gissar  and  Zeravshan  Ranges),  Tadjik  Depression  and  western  Pamir.  This  is 
the  most  active  part  of  Soviet  Central  Asia  and  has  been  extensively  studied  by  Soviet  groups  at  the 
Tadjik  Institute  of  Seismoresistant  Construction  and  Seismology  (TISSS)  in  Dushanbe  and  the 
Complex  Seismological  Expedition  (KSE)  in  Garm,  an  outpost  of  the  Insititute  of  the  Physics  of  the 
Earth,  Moscow.  It  is  also  the  site  of  our  earlier  work  on  induced  seismicity  at  Nurek  Reservoir 
(Simpson  and  Negmatullaev  81,  Keith  et  al  82,  Leith  et  al  81)  and  regional  tectonics  (Kristy  and 
Simpson  1980,  Leith  and  Alvarez  1985,  Leith  and  Simpson  1986a,b).  Figure  15  identifies  in  more 
detail  some  of  the  geographic  features  in  a  subsection  of  Figure  14  near  Nurek  reservoir. 

Figure  16  represents  topography  with  intensity  directly  proportional  to  elevation.  In  Figure  17  a 
filter  has  been  used  to  locally  enhance  the  contrast  within  different  topographic  provinces.  Thus  the 
relationship  between  intensity  and  elevation  is  no  longer  linear  or  spatialy  constant,  but  many  of  the 
subtle  topographic  features  are  more  clearly  distinguished.  In  Figure  18  a  simple  shaded  relief  map 
has  been  produced  by  taking  the  first  spatial  derivative  of  the  elevation  in  the  NW-SE  direction,  giving 
the  impression  of  illumination  from  the  NW.  This  filter  accentuates  NE-SW  striking  features  (eg. 
ridges  in  the  Tadjik  Depression).  A  similar  image,  with  illumination  from  the  NE  (Figure  19) 
accentuates  features  perpendicular  to  those  in  Figure  18. 

One  of  the  first  order  features  obvious  in  Figures  16- 19  is  the  contrast  between  the  topographic 
style  in  the  northern  half  of  the  image  and  that  in  the  south.  The  northern  segment,  with  higher  and 
smoother  average  topography,  a  general  EW  fabric  and  deeply  incised  rivers  is  composed  primarily  of 
Paleozoic  granites  of  the  Southern  Tien  Shan.  The  southern  part,  with  pronounced,  sharp  ridges 
oriented  more  NE-SW,  shows  the  results  of  recent  folding  and  thrusting  of  the  young  foreland  fold 
and  thrust  belt  of  the  Tadjik  Depression  (Leith  and  Alvarez,  1985).  The  seismicity  of  this  area  (Figure 
20)  and  its  relationship  to  the  major  physiographic  provinces  is  described  in  detail  by  Leith  and 
Simpson  (1986a).  Three  clear  types  of  seismicity  are  found  in  the  shallow  earthquakes  of  the  region 
(in  addition  to  the  intermediate  depth  activity  of  the  Hindu  Kush).  The  infrequent,  large  earthquakes 
(up  to  M=8+)  are  confined  to  the  Gissar  Kokshal  zone,  along  the  front  of  the  Southern  Tien  Shan,  but 
that  zone  has  a  relatively  low  level  of  microearthquake  activity.  The  greatest  concentration  of 
seismicity,  in  terms  of  numbers,  is  within  the  sediments  of  the  Tadjik  Depression,  but  there  are  few 
events  here  with  magnitude  greater  than  5.5.  Beneath  the  sediments,  the  third  type  of  seismicity  is  that 
within  the  basement  itself.  The  activity  in  the  Depression  is  especially  high  in  its  northeastern  end  near 


Garm,  in  the  area  of  greatest  convergence  between  the  Pamir  and  Tien  Shan.  In  this  area,  the 
sediments  have  been  thrust  up  and  out  of  the  Depression  itself,  onto  the  Gissar  Kokshal  zone, 
producing  a  complex  three  dimensional  distribution  of  seismicity  in  which  all  three  types  of  activity 
overlap  (Leith  and  Simpson,  1986a). 
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A  strong  concentration  of  shallow  seismicity  (up  to  magnitude  5)  stands  out  in  the  relatively  low 
activity  of  the  southern  Tadjik  Depression  near  the  town  of  Kulyab  (Figure  20).  There  are  also  large 
salt  domes  exposed  at  the  surface  in  this  same  region.  We  have  used  Landsat  and  topographic  data  in 
conjunction  with  the  seismicity  catalog  to  investigate  the  correlation  between  the  seismicity  and  the  salt 
domes  (Leith  and  Simpson,  1986b). 

Figure  21  shows  a  black  and  white  version  of  the  Landsat  image  for  the  Kulyab  area.  Areas  of 
high  infrared  return  (vegetation)  have  been  set  black  in  this  image  so  that  the  dark  regions  in  valleys 
are  irrigated  farm  land.  The  two  light  colored  circular  features  are  salt  domes.  Figure  22  shows  a 
perspective  view  (from  the  SW)  derived  by  creating  a  3D  image  from  the  digital  terrain  data  and  adding 
texture  from  the  Landsat  image.  The  major  dome  to  the  south  (Mumin  dome)  is  devoid  of  vegetation 
and  stands  870  meters  above  the  surrounding  plain.  The  high  topography  beyond  the  Mumin  dome  is 
associated  with  the  Sari-Chashma  salt  diapir,  a  buried  structure  identified  in  geophysical  surveys. 
Figure  23  shows  the  original  and  shaded  relief  versions  of  the  topographic  data  with  seismicity 
superimposed.  Most  of  the  cluster  of  seismicity  in  the  southeastern  comer  of  the  image  followed  a 
m=5.2  earthquake  on  April  2,  1973. 

It  is  our  conclusion  (Leith  and  Simpson,  1986b)  that  the  seismicity  is  clearly  related  to  the  salt 
doming.  The  highest  level  of  current  activity  is  not  located  at  one  of  the  domes  which  has  pierced  the 
surface,  but  is  near  the  buried  Sari-Chashma  diapir,  a  buried  structure  which  appears  to  be  actively 
rising.  Because  of  the  lack  of  depth  control  in  the  Soviet  hypocentral  locations,  it  is  not  clear  whether 
the  earthquakes  result  from  increased  stress  within  the  sediments  above  the  dome  structure  or  are 
related  to  fracturing  near  or  within  the  salt  layer  at  depth.  To  our  knowledge,  this  is  the  only  reported 
case  of  seismicity  associated  with  active  salt  doming.  More  detail  is  to  be  found  in  Leith  and  Simspon 
(1986b,  included  as  appendix  to  this  report). 
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Figure  3  Cumulative  and  annual  numbers  of  earthquakes  in  the  ZSSSR  and  ISC  catalogs. 


Figure  6  Difference  in  latitude  and  longitude  (dy  and  dx)  between  ZSSSR  location  (origin  in  each 
graph)  and  ISC  location  for  all  earthquakes  common  to  both  catalogs,  as  a  function  of  the  number  of 
stations  used  for  location  by  ISC. 


gure  8  Difference  in  depth  (dz)  between  ISC  and  ZSSSR  catalogs  as  a  function  of  the  number 
stations  used  for  location  by  ISC. 


Figure  9  Difference  in  origin  time  (dt)  between  ISC  and  ZSSSR  catalogs  as  a  function  of  the 
number  of  stations  used  for  location  by  ISC. 


DEEP  EVENTS 


Figure  1 0  Difference  in  origin  time  (dt)  between  ISC  and  ZSSSR  catalogs,  as  a  function  of  the 
number  of  stations  used  for  location  by  ISC,  for  deep  events  (those  for  which  both  the  ZSSSR  and 
ISC  depths  are  greater  than  60  km).  (Note  -  positive  dt  is  for  ISC  origin  times  later  than  ZSSSR). 


SHALLOW  EVENTS 

(ZOSCXIOO  ;  ZIS0VX60) 


is  Figure  10,  but  for  shallow  events  (those  for  which  the  ZSSSR  depth  is  less 
ISC  depth  is  less  than  60  km,  but  neither  are  the  default  depths  of  0  or  33  km). 


Figure  12  Difference  in  origin  time  (dt)  vs  difference  in  depth  (dz)  for  shallow  events  (as  defined  in 
Figure  11). 


3.29 


Figure  1 3  ZSSSR  energy  class  (k)  vs  ISC  magnitude  (m^)  and  number  of  earthquakes  in  ISC  for 

each  0. 1  magnitude  increment  Least  squares  regression  (with  weight  inversely  proportional  to  square 
of  the  standard  deviation  of  the  average  magnitude  for  each  magnitude  interval)  is  shown  and 
compared  to  the  standard  formula  of  Rautian. 
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14  (from  Simpson  and  NegmatuUaev,  1981). 


Figure  19  Shaded  relief  map  with  illumination  from  the  northeast.  Comparison  with  Figures  14  and 
1 8,  especially  of  river  valleys  in  the  northeast  segment,  shows  how  different  gradient  determinations 
(illumination  directions)  enhance  features  at  different  azimuths. 


Figure  22  Perspective  view  of  the  area  show  in  Figure  21.  Digital  terrain  data  have  been  used  to 
create  a  pespective  image  (viewed  from  the  southwest)  and  the  Landsat  image  has  been  registered  with 
it  to  create  surface  texturing.  The  salt  dome  in  the  lower  right  stands  more  than  800  m  above  the 
surrounding  valley. 
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Figure  23  (top)  Raw  digital  terrain  data  for  the  Kulyab  area,  (middle)  Shaded  relief  map  with 
seismicity  superimposed,  (bottom)  Location  map  (from  Leith  and  Simpson,  1986). 
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EARTHQUAKES  RELATED  TO  ACTIVE  SALT  DOMING 
NEAR  KULYAB,  TADJIKISTAN,  USSR 

William  Leith  and  David  W.  Simpson 
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Abstract.  In  the  region  near  Kulyab,  Tadjikistan, 
hundreds  of  shallow  earthquakes  with  magnitudes 
greater  than  234  have  been  reported  in  Soviet  yearly 
catalogs  since  1961.  This  area  appears  as  a  well- 
defined  cluster  of  activity,  distant  from  the  line  of 
epicenters  that  defines  the  Gissar-Kokshal  seismic 
zone,  to  the  north  of  the  Pamir  ranges.  The  geology 
of  this  region  is  dominated  by  the  presence  of 
numerous  salt  domes,  surrounded  by  Neogene  and 
Quaternary  continental  deposits.  The  spatial  relation¬ 
ship  between  these  earthquakes  and  the  salt  domes 
suggests  that  the  two  phenomena  may  be  related. 
Moderate  earthquakes  (M  >  5)  occurred  in  1972  and 
10711.  and  intensities  of  surface  shaking  greater  than 
MM=6  were  reported  from  earthquakes  in  19.37,  1952, 
1969.  1972,  1973  and  1978.  The  earthquake  on  2  April 
1973  and  its  aftershocks  were  located  in  a  region 
where  no  salt  domes  have  been  mapped  at  the  sur¬ 
face.  However,  a  buried  salt  diapir  has  been  mapped 
at  depth  by  geophysical  means.  These  earthquakes 
may  result  from  active  salt  diapirism  at  depth.  The 
mechanism  for  producing  this  seismicity  could  be 
either  the  active  fracturing  of  the  cap  rock  by  the  ris¬ 
ing  diapir,  or  the  concentration  of  tectonic  stresses  in 
the  thinned  section  above  and  adjacent  to  the  diapir. 
The  salt-related  earthquakes  may  produce  lower  fre¬ 
quency  radiation  than  other  events  of  the  same  size. 

Introduction 

Salt  deposits,  and  specifically  salt  domes,  are 
currently  of  special  interest  as  possible  sites  for 
decoupling  of  underground  nuclear  explosions  [Han¬ 
non,  1985],  or  for  the  underground  storage  of  radioac¬ 
tive  nuclear  waste,  petroleum  or  gas  condensate 
[Cohen,  1977;  Borg,  1983|.  The  relatively  structure¬ 
less,  homogeneous  nature  of  salt  formations,  com¬ 
bined  with  the  high  solubility  of  rock  salt,  lend  these 
formations  to  the  mining  of  the  large  stable  cavities 
that  would  be  required  to  decouple  a  small  under¬ 
ground  nuclear  explosion.  In  contrast  with  bedded 
salt  deposits,  salt  domes  are  considered  ideal  cavity 
sites  because  of  their  large  vertical  dimensions  and 
relatively  pure  salt  content.  In  terms  of  explosion 
decoupling  and  the  detection  of  clandestine  nuclear 
weapons  tests,  earthquakes  that  are  spatially  associ¬ 
ated  with  salt  domes  pose  the  obvious  problem  of 
discrimination.  In  terms  of  engineering  activities, 
either  on  the  surface  or  within  the  salt  dome,  earth¬ 
quakes  represent  a  hazard  to  the  stability  of  any  rav- 
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ity  constructed.  Also,  as  potential  sites  for  hazardous 
waste  or  other  liquid  storage,  they  may  affect  the 
short-term  permeability  of  the  cavity  walls. 

There  is  an  extensive  literature  on  salt  diapirism 
[e.g.,  Lerche,  1986).  Salt  domes  are  produced  by  the 
post-depositional  flow  of  bedded  salt  upward  through 
the  overlying  rock  formations,  because  of  the  gravity 
(density)  instability  inherent  in  buried  salt.  It  is  gen¬ 
erally  considered  necessary  for  salt  dome  regions  to 
have  undergone  some  degree  of  tectonism  (associated 
with  sedimentary  loading)  in  order  for  diapirs  to 
develop.  Therefore,  the  process  of  salt  dome  forma¬ 
tion  is  unstable  and  episodic  and  the  rates  of  salt 
dome  uplift  are  highly  variable,  depending  on  both 
the  tectonics  of  the  region  and  the  state  of  maturity 
of  the  diapir  [Jackson  and  Talbot,  1986]. 

Salt  is  an  extremely  weak  rock,  and  deforms  duc- 
tilely  at  geological  strain  rates  [Carter  and  Hansen, 

1 083].  Salt  areas  are  generally  not  associated  with 
earthquakes,  and  to  date  there  is  no  literature  associ¬ 
ating  earthquakes  with  active  salt  doming.  Where 
studies  have  been  made  of  the  seismicity  of  salt  dome 
areas  in  East  Texas  and  the  Gulf  Coast  (Dorman  and 
others,  1977;  Racine  and  Klouda,  1979],  no  definitive 
correlations  between  salt  structures  and  earthquakes 
have  been  made.  Nevertheless,  sounds  of  salt  move¬ 
ment  have  been  detected  in  salt  dome  regions  in  Iran 
[Kent,  1979],  and  rock  bursts  have  occurred  in  under¬ 
ground  salt  mines  [see  Baar,  1977],  Those  salt  dome 
areas  where  detailed  seismicity  studies  have  been 
made  (east  Texas,  Louisiana,  Oklahoma)  are  also 
areas  where  the  domes  are  mature  [Jackson  and  Tal¬ 
bot,  1986)  and  tectonic  strain  rates  are  low.  Our 
study,  from  an  area  that  is  currently  deforming  at  a 
relatively  high  rate,  therefore  contrasts  with  studies 
of  the  seismicity  of  salt  dome  areas  in  the  U.S. 

This  letter  describes  the  geological  and  seismologi- 
cal  attributes  of  the  Kulyab  salt  dome  area,  in  the 
southwestern  Tadjik  SSR.  The  seismicity  that  marks 
this  area  stands  out  as  a  cluster  of  epicenters  distinct 
from  the  main  trend  of  the  Gissar-Kokshal  seismic 
zone  (figure  1),  50-100  km  north  of  Kulyab  [Leith 
and  Simpson,  1986). 

Geology 

The  Soviet  republic  of  Tadjikistan  is  marked  by  an 
active  tectonics  and  a  high  rate  of  natural  seismicity. 
As  part  of  the  India-Asia  collision,  the  Jurassic  and 
younger  sediments  that  fill  the  sedimentary  basin  of 
the  Tadjik  Depression,  a  relatively  low  region  that 
forms  most  of  southwest  Tadjikistan  (Figure  2),  have 
been  deformed  by  folding  and  thrusting  in  the  zone 
between  the  I’amir  and  Tieti  Shan  ranges. 
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Fig.  1.  Shallow  earthquakes  for  the  area  of 
southwestern  Tadjikistan  from  1970-1980.  Seismicity 
is  most  intense  within  the  frontal  (northern)  portions 
of  the  Vakhsh  fold/thrust  belt  and  along  a  portion  of 
the  Darvaz  fault.  The  cluster  of  seismicity  in  the 
lower  portion  of  the  figure  (arrow)  occurs  near  the 
salt  domes  mapped  in  Fig.  2,  and  well  away  from  the 
belt  of  epicenters  that  marks  the  fold/thrust  belt.  A 
detail  of  this  area  is  shown  in  Fig.  3. 


Throughout  the  Depression,  Jurassic  evaporites  form 
the  base  of  the  sections  exposed  along  thrust  faults  at 
the  surface,  and  the  overall  structure  of  the  Depres¬ 
sion  suggests  that  these  basal  evaporites  form  the 
decollement  across  which  the  thrusts  are  displaced 
[Leith  and  Alvarez,  1986],  Hasin  formation  in  mid- 
Jurassic  time  was  followed  by  the  deposition  of  large 
thicknesses  of  evaporites.  Since  the  end  of  the  Juras¬ 
sic,  shallow-water  marine  sedimentation  in  the 
Depression  proceeded  at  a  relatively  slow  rate,  as  the 
basin  subsided  following  the  rifting  event  (Leith, 
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Fig.  2.  Generalized  crustal  structure  and  geology  of 
the  northern  Tadjik  Depression.  Salt  domes  are 
located  near  Kulyab  in  the  southeastern  Tadjik 
Depression,  10-50  km  from  the  Soviet- Afghan  border, 
which  follows  the  Pyand/h  river  (dash-dot  line). 
Depth  to  basement  beneath  the  Mesozoic  and  Oeno 
zoic  sediments  of  the  Tadjik  Depression  is  contoured 
in  3  km  intervals  (thin  lines). 
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Fig.  3.  (a)  Earthquake  epicenters  from  the  Soviet 
Catalogs  for  the  Kulyab  region,  mapped  on  shaded 
digital  terrain  data,  (b)  Sketch  shows  salt  deposits, 
major  faults  and  geography.  There  is  a  close  spatial 
correlation  between  some  groups  of  earthquakes'  and 
mapped  salt  domes.  The  cluster  of  earthquakes  in  the 
lower  right  is  associated  with  the  buried  Sari 
Chashma  salt  diapir  (detail  in  Fig.  1).  Location  accu¬ 
racy  for  most  events  is  listed  as  ±  10  km. 


1985].  Adjacent  to  the  Pamir  range  in  the  region  of 
Kulyab.  this  was  followed  in  about  middle  Oligocene 
time  by  the  deposition  of  more  than  a  km  of  Neogene 
molasse  in  a  north-south  trending  trough  of  thick 
and,  now,  gently-folded  rocks  (the  Kulyab  synclino- 
riuin).  This  rapid  loading  of  the  Jurassic-Paleogcne 
section  has  provided  the  drive  for  the  salt  migration. 

The  Kulyab  area,  along  the  western  liinh  of  the 
Kulyab  synclinorium.  is  marked  by  the  outcrop  of 
several  large  bodies  of  Jurassic-Lower  Cretaceous  salt 
[Luchnikov,  1982].  Some  are  obvious  domes  (Figure 
3),  while  others  mark  the  traces  of  mapped  thrust 
faults.  Several  of  the  domes  have  prominent  topo¬ 
graphic  relief.  For  example,  at  the  Mumin  dome, 
south  of  Kulyab  (see  Figures  3,  -1),  the  central  salt 
stock  rises  more  than  870  m  above  the  surrounding 
plain.  This  feeds  an  oblate  glacial  salt  flow  above 
folded  and  faulted  upper  Tertiary  strata  [Sadikov, 
1982j. 

Geophysical  studies  and  exploratory  drilling  that 
have  accompanied  the  search  for  oil  and  gas  in  the 
Kulyab  area  have  helped  to  delineate  the  distribution 
of  salt  at  depth  [Luchnikov,  1982:  Azimov  and  others, 
1982|,  and  have  enabled  the  mapping  of  a  zone  of 
overpresMired  formations  in  the  west  of  Kulyab 
[Kalomazov  and  Vakhlikov,  I975j.  A  well  drilled  in 
the  northern  dome  in  Figure  3  (‘‘Sartis")  encountered 
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Fig.  4.  Epicenters  from  the  April  2nd,  1973  earth¬ 
quake  and  its  aftershocks,  mapped  on  a  Landsat 
image  of  the  Kulyab  area.  The  location  of  a  mapped 
subsurface  diapir  (dashed  line)  is  from  Kon’kov 
(1974],  This  earthquake  sequence  apparently 
occurred  at  shallow  depth  along  the  northern  margin 
of  the  diapir. 

pressures  of  1.6  and  1.8  times  hydrostatic  in  Turanian 
age  strata  at  depths  of  2642  and  2695  m,  respectively. 
In  general,  the  westernmost  (up-structure)  Helds  are 
overpressured  (up  to  twice  hydrostatic)  while  those  to 
the  east  are  not. 

Seismicity 

Epicenters  of  earthquakes  in  the  Kulyab  area  from 
1964-1980,  from  the  yearly  catalogs,  Zemletracenniye 
v  SSSR  (hereafter,  ‘ZSSSR  catalogs’),  are  shown  in 
Figure  .'I,  along  with  shaded  digital  topography. 
There  are  352  earthquakes  listed  in  the  catalogs 
(many  duplicate  locations  do  not  appear  in  the  fig¬ 
ure),  all  with  magnitudes  greater  than  about  2.5.  All 
earthquakes  in  this  area  are  shallow  (<  20  km),  and 
there  is  a  clear  association  of  some  events  with  the 
mapped  salt  domes.  The  large  cluster  of  activity  in 
the  lower  right  corner  of  the  figure  occurred  following 
the  M=5.2  earthquake  on  02  April  1973  (hereafter, 
“1973  earthquake").  This  sequence  is  associated  with 
a  diapir  at  depth,  discussed  IipIow. 

(A  Soviet  regional  catalog  lists  some  smaller  events 
for  the  Kulyab  area  (to  M  <  2)  but  is  incomplete.  A 
map  showing  some  smaller  events  was  published  in 
the  1979  ZSSSR  ealalog  (p.  30  overleaf),  but  we  can¬ 
not  specify  location  accuracy  and  salt  structures  are 
not  located  on  the  map). 

For  the  region  of  Figure  3,  the  b-value  is  approxi¬ 
mately  1. 1-1.2  for  earthquakes  from  1964-1980.  This 
is  relatively  high  in  comparison  with  the  western  por¬ 


tion  of  the  Vakhsh  fold/thrust  belt,  and  is,  in  fact, 
comparable  to  that  of  the  reservoir-induced  seismicity 
at  Nurek  [Simpson  and  Negmatullaev,  1981).  Induced 
earthquakes  are  triggered  by  effective  stress  changes 
resulting  from  changes  in  the  subsurface  pore  pres¬ 
sure  regime  [Simpson,  1986],  and  relatively  high  b- 
vaiues  are  typical  of  induced  seismic  activity  [Gupta 
and  Rastogi,  1976).  We  suggest  that  the  high  b-value 
that  we  have  determined  for  the  Kulyab  region  may 
be  “normal”  for  that  area,  because  of  the  existence  of 
the  naturally  high  formation  pressures  noted  above. 
The  high  b-values  may  thus  be  a  fluid  pressure  effect 
that  distinguishes  this  region  from  other  areas  in 
western  Tadjikistan. 

The  Kulyab  area  has  produced  several  earthquakes 
with  magnitudes  of  about  5  in  this  century  [Gubin, 
1960),  including  those  in  1937,  1952,  1959,  1972,  1973 
and  1978.  From  these,  the  April  2nd,  1973  earthquake 
(M«=k5.2)  is  the  best  studied  to  date.  Kon’kov  [1974] 
combined  data  from  the  Kulyab  seismic  station  with 
til Hcroftcislil }r  data  to  relocate  the  epicenter  of  the 
main  shock,  which  lies  15  km  NNW  of  the  epicenter 
reported  in  ZSSSR,  at  a  depth  of  12  km.  It  lies  near 
the  northern  edge  of  the  buried  Sari-Chashma  diapir, 
10  km  SE  of  Kulyab.  The  sequence  is  plotted  with 
geologic  and  geographic  data  in  Figure  4.  The  earth¬ 
quakes  apparently  occurred  along  the  northern 
margin  of  the  Sari-Chashma  diapir,  at  depths  from  0 
to  10  km.  The  main  event  (April  2,  00h02m43s)  was 
included  in  a  spectral  study  of  Central  Asian  earth¬ 
quakes  by  Zapol’skiy  and  Loginova  (1984;  figure  2), 
which  suggests  that  this  event  produced  lower  fre¬ 
quency  seismic  waves  than  other  earthquakes  of  com¬ 
parable  size.  Low  frequency  microearthquakes  near 
Nurek  reservoir  may  also  be  associated  with  salt 
deposits,  but  more  work  needs  to  be  done  to  deter¬ 
mine  the  cause  of  the  spectral  anomaly. 

Discussion 

Salt  diapirism  in  active  tectonic  regions  can  be 
associated  with  very  high  rates  of  deformation.  It  is 
therefore  more  appropriate  to  compare  the  dynamics 
of  the  Kulyat  salt  domes  with  those  of  the  active  salt 
dome  area  of  southern  Iran  [Ala,  1974],  rather  than 
the  more  stable  region  of  extensive  salt  domes  along 
the  U.S.  coast  of  the  Gulf  of  Mexico.  For  example,  a 
detailed  study  of  salt  dynamics  has  been  made  for  the 
Kuh-c-Namak  intrusion  in  Iran  [Talbot  and  Jarvis, 
1984].  This  salt  stock  maintains  the  highest  known 
natural  rate  of  salt  extiusion,  estimated  by  Talbot 
and  Jarvis  [1984]  to  average  170  mm/yr  (or  about  1 
km/6000  yr).  This  represents  an  extremely  high 
strain  rate  in  the  salt  stock,  and  a  significant  addi¬ 
tional  shear  stress  in  the  adjacent  country  rock  [D. 
Davis,  pers.  comm.,  1986). 

Because  of  the  very  artive  nature  of  the  tectonics  of 
Tadjikistan,  we  suspect  that  the  diapirs  near  Kulyab 
reflect  the  same  high  degree  of  activity  as  the  Iranian 
intrusions.  For  comparison,  the  Mumin  dome  (which 
is  more  than  7  km  wide  and  rises  more  than  870  m 
above  the  adjacent  alluvial  plain;  see  Fig.  4)  is 
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approximately  the  same  size  as  the  Kuh-e-Namak 
intrusion,  and  the  domes  are  both  in  similar  arid 
environments.  The  Murnin  dome  must,  therefore, 
maintain  a  rate  of  salt  extrusion  of  the  order  of  the 
Kuh-e-Naniak,  in  order  to  maintain  comparable 
topography.  Unfortunately,  data  on  the  deep  struc¬ 
ture  of  the  Sari  Chashma  diapir  or  surface  strains  in 
the  local  area  of  the  diapir  are  unavailable.  Thus,  we 
are  currently  not  able  to  determine  the  amount  or 
character  of  the  deformation  that  is  associated  with 
the  growth  of  this  salt  body. 

Fdsewhere  in  the  Tadjik  Depression,  shallow  seismi¬ 
city  is  largely  confined  to  the  sedimentary  section 
above  the  salt-layer  decollement  (Leith  and  Simpson, 
1986],  and  salt  is  generally  thought 'too  weak  to  hold 
the  large  elastic  strains  that  would  be  necessary  to 
produce  a  moderate-size  earthquake.  Thus,  we 
suspect  that  the  1974  earthquake  sequence  reflects 
fracturing  of  the  country  rock  above  or  adjacent  to 
the  Sari  Chashma  diapir,  and  not  the  seismic  release 
of  stored  elastic  strains  within  the  salt  itself.  How¬ 
ever,  depths  and  focal  mechanisms  of  earthquakes  in 
the  Kulyab  area  are  poorly  constrained,  and  we  must 
await  further  data  before  speculating  on  the  mechan¬ 
ics  of  the  deformation. 
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1.  Introduction 


Apart  from  their  utility  to  the  study  of  regional  seismicity,  the  use  of  regional  seismic 
networks  for  the  study  of  the  properties  of  teleseismic  travel  times  and  waveforms  is  well 
established.  Studies  of  the  crust  and  upper  mantle  subjacent  to  array  installations  have  been  done  at 
LASA  and  NORSAR  [e.g.  Aki,  1973;  Aki  et  al.,  1977]  and  SCARLET  [e.g.  Humphreys  et  al., 
1984].  Arrays  have  also  been  useful  for  the  study  of  distal  structure.  Bungum  and  Capon  [1974] 
used  LASA  recordings  to  measure  apparent  back  azimuths  of  Rayleigh  waves.  Engdahl  [1968] 
used  LASA  records  to  identify  and  study  PKiKP  and  PmKP  arrivals.  Husebye  et  al.  [1976]  used 
NORSAR  to  measure  the  properties  of  PKIKP  precursors,  and  Chang  and  Cleary  [1978,  1981] 
and  Doornbos  [1980]  used  arrays  to  measure  the  properties  of  PKKP.  NORSAR  has  also  been 
used  to  study  the  nature  of  P-diffracted  and  its  decay  into  the  core  shadow  [ Doornbos  and  Mondt, 
1979],  The  recently  operational  NARS  network  [ Dost  et  al.,  1984]  is  being  used  to  assemble  a 
broad  range  of  observations  of  fundamental  and  higher-mode  surface  waveform  characteristics  and 
to  study  the  properties  of  phases  converted  at  upper-mantle  discontinuities  [Paulssen,  1985;  Dost, 
1986].  In  this  report,  we  explore  the  utility  of  the  New  York  State  Seismic  Array  (NYSSA), 
operated  by  Lamont-Doherty,  for  studies  of  regional  variability  of  structure  and  the  properties  of 
the  core-mantle  boundary  and  D". 

Waveform  studies  using  arrays  offer  the  advantage  of  signal  enhancement  by  assuming  (with 
some  justification)  that  the  background  noise  is  uncorrelated  at  frequencies  of  interest.  Various 
stacking  and  beaming  procedures  have  been  designed  to  optimize  the  response  of  the  array  to  the 
particular  signal  being  evaluated  [e.g.  Capon  et  al.,  1967;  Nolet  and  Panza,  1976;  Husebye  et  al., 
1976;  Ingate  et  al.,  1985].  Such  procedures  depend  critically  upon  the  physical  design  of  the  array, 
that  is,  the  disposition  of  the  instruments.  While  these  procedures  elevate  the  accuracy  of  various 
sorts  of  signal  measurements,  of  somewhat  more  importance  for  global  studies  is  the  position  of 
the  array  relative  to  seismogenic  regions  [Gee  et  al.,  1985].  NYSSA  is  situated  ideally  with  respect 
to  the  seismogenic  zones  of  the  southwestern  Pacific  for  the  study  of  the  core  phases  belonging  to 
the  PKP  group.  In  addition  to  their  instrinsic  interest,  PKP  phases  arrive  at  the  array  with  nearly 
vertical  incidence,  and  provide  a  reasonable  input  signal  with  which  near- receiver  transfer  function 
can  be  evaluated,  and  small-scale  regional  structure  can  be  estimated.  Of  importance  to  the  accurate 
interpretation  of  these  arrivals  is  a  measure  of  the  coherence  of  the  signal,  and  the  response  of  the 
array  to  incoming  signals  of  near  vertical  incidence.  In  this  report,  we  assess  the  utility  of  using 
NYSSA  for  the  analysis  of  these  phases.  Initial  estimates  of  spatial  coherence  over  a  restricted 
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Figure  1:  Configuration  of  New  York  State  Seismic  Array  as  of  December  10,  1986. 
LDGO  is  located  near  station  PAL  in  southeastern  New  York  State.  Roman  numerals 
refer  to  subnetworks;  each  will  have  a  central  data  recording  /  transmission  site  tied  to 
LDGO. 
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portion  of  the  array  will  be  presented.  These  measurements  are  important  to  a  thorough 
understanding  of  higher-frequency  regional  wave  propagation,  a  program  of  interest  to  DARPA. 
Future  DARPA-funded  research  at  LDGO  will  be  concerned  with  the  robust  estimation  of  spatial 
coherence  using  NYSSA  and  specially  installed  sub-arrays,  and  the  implications  for  crustal  and 
mande  heterogneity  beneath  the  array. 

2.  Network  configuration 

The  New  York  State  Seismic  Array  (NYSSA)  was  developed  and  installed  by 
Lamont- Doherty  Geological  Observatory  (LDGO)  to  monitor  seismicity  in  the  northeastern  United 
States.  The  current  configuration,  shown  in  Figure  1,  comprises  27  stations  partitioned  into  four 
sub-regional  networks  in  northern  New  York  State,  the  Adirondacks,  southeastern  New  York  and 
northern  New  Jersey,  and  western  New  York.  Station  sites  were  chosen  to  provide  optimal  array 
characteristics  for  local  and  regional  seismic  monitoring  experiments.  The  dotted  line  surrounds  the 
region  in  which  there  is  nearly  complete  detection  of  events  above  magnitude  2.0.  Sensors  have 
fundamental  periods  of  one  or  two  seconds;  together  with  associated  electronics,  the  total  system 
response  is  flat  to  acceleration  between  .1  and  1  Hz  and  flat  to  velocity  between  1  and  10  Hz 
(Figure  2).  System  output  is  digitized  using  a  12-bit  A/D  with  100  Hz  sampling,  giving  a  usable 
response  to  30  Hz.  The  pre-event  memory  is  variable,  and  can  be  adjusted  to  a  maximum  of  68 
seconds. 


NYSSA  Nominal  Response  to  Displacement 


Figure  2:  Nominal  NYSSA  instrument  displacement  response,  alter  Lee  and  Stewart 
[1981], 


The  network  operates  in  triggered  mode  using  a  short-term  average  -  long-term  average 
algorithm  tuned  to  seismicity  levels  in  the  Northeast.  Detected  events  are  transmitted  over  leased 


telephone  lines  to  LDGO  using  analog  telemetry.  Although  the  triggering  algorithm  is  designed  to 
detect  local  and  regional  events  reliably,  teleseismic  arrivals,  particularly  those  having  higher 
frequency  content,  will  trigger  the  network  and  be  recorded.  Since  the  network  trigger  must  be 
tuned  to  avoid  saturation  on  small  local  events,  the  effective  magnitude  threshold  for  teleseisms  at 
core-study  distances  is  greater  than  about  mb  =  6.0.  Selected  analog  recordings  have  been  made  of 
events  with  body-wave  magnitudes  as  low  as  5.2  -  5.5,  however,  with  sufficient  "eyeball" 
signal-to-noise  levels  to  be  useful  for  core-phase  analysis.  To  trigger  on  these  events,  and  to  avoid 
saturation  by  small  local  events,  the  triggering  algorithm  must  be  modified  to  make  near  real-time 
estimates  of  incoming  signal  azimuth  and  incidence  angle. 

3.  Comparison  of  NYSSA  and  NORSAR  for  PKP  analysis 


The  New  York  State  Seismic  Array  (NYSSA)  was  established  principally  to  monitor  regional 
seismicity  in  the  northeastern  United  States.  Teleseismic  triggers  typically  have  been  used  to 
provide  external  checks  on  instrument  polarities.  NORSAR,  on  the  other  hand,  was  established 
principally  as  a  network  with  a  well-defined  array  response  to  incoming  teleseismic  arrivals  in  order 
to  resolve  questions  associated  with  nuclear  event  detection  and  discrimination,  though  it  has 
proven  very  successful  in  addressing  basic  earth  structure  problems  [e.g.  Husebye  et  al.,  1976;  Aki 
etal.,  1977 ;  Door nbos  and  Mondt,  1979]. 


NYSSA  and  NORSAR  Deltas  for  1985 


Figure  3:  Histogram  of  expected  arrivals  per  sensor  at  NYSSA  (squares)  and  NORSAR 
(triangles)  for  Pacific  Rim  earthquakes  with  mb  >  5.5.  A  5°  running  mean  has  been 
applied  to  simulate  wider  aperture.  Seismicity  is  taken  from  1985  PDE  catlaogs. 


Figure  3  shows  a  comparison  of  expected  annual  arrivals  (per  sensor)  at  NYSSA  and 
NORSAR  as  a  function  of  epicentral  distance  from  Pacific  rim  events  with  mh  greater  than  5.5, 
using  events  in  the  1985  PDE  catalog.  The  number  of  arrivals  integrated  over  the  distance  range 
1 10°  -  160°  is  about  the  same  for  both  arrays,  but  the  distributions  are  different.  The  NYSSA  array 
has  a  peak  in  the  arrival  histogram  between  1 17°  and  about  130°,  with  roughly  double  the  number 
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of  expected  arrivals  at  NORSAR,  while  NORSAR  should  see  about  50%  more  arrivals  bewteen 
130°  and  140°  than  NYSSA.  (Indeed,  the  study  of  DF-precursors  by  Husebye  et  al.  11976]  used 
events  concentrated  around  131°  and  136°,  corresponding  roughly  to  the  peaks  in  the  NORSAR 
histogram  in  Figure  3.)  Overall  the  NYSSA  histogram  is  smoother,  reflecting  a  more  favorable 
distribution  of  event  geometry  and  a  wider  aperture.  NORSAR,  however,  should  have  greater 
detection  capabilities  due  to  its  more  controlled  array  response  [Dahle  et  al.,  1975],  a  factor  which 
we  will  investigate  for  NYSSA.  The  total  azimuthal  sampling  is  comparable  between  the  two 
arrays,  though  of  course,  the  patches  of  the  core  sampled  by  each  array  are  different.  Our  point 
here  is  not  to  focus  on  dissimilar  aspects  of  NYSSA  and  NORSAR,  but  merely  to  demonstrate  that 
both  NORSAR  and  NYSSA  are  potentially  capable  of  illuminating  equivalent  portions  of  P' 
branches.  Examples  of  DF-precursor  waveforms  from  NORSAR  are  published  in  Husebye  et  al. 
[1976].  Examples  of  waveforms  from  NYSSA  are  shown  in  the  next  section. 

4.  Initial  observations  of  PKP 

We  have  examined  22  teleseismic  triggers  from  events  in  the  southwest  Pacific  and  Indonesia 
covering  an  epicentral  range  from  1 1 1°  to  148°.  The  network  triggers  on  P'DF  or  P'Bcfor  these 
events,  and  there  is  enough  pre-event  memory  in  the  standard  configuration  to  obtain  good 
recordings  of  DF-precursors  or  the  AB  -  BC  -  DF  crossovers.  At  these  azimuths,  the  network 
aperture  is  about  5°.  Examples  of  three  record  sections  exhibiting  DF-precursors  are  shown  in 
Figures  4,  5,  and  6.  Figure  7  shows  a  record  section  illuminating  the  BC  -  DF  crossover. 

Event  860115.2017,  shown  in  Figure  4,  is  one  of  the  closest  events  for  which  we  have 
observed  unambiguous  precursors  (1 19°  -  124°).  (To  our  knowledge,  these  are  among  the  closest 
observations  of  DF-precursors  yet  collected.)  DF-precursors  are  well-observed  on  stations  beyond 
122°  (these  are  Adirondack  sites)  with  impulsive  onset  and  a  ringing  behavior  of  nearly  constant 
amplitude  up  to  the  DF  arrival.  Total  length  of  the  precursive  wavetrain  is  1 1  -  12  seconds.  Noise 
levels  at  the  stations  vary,  and  the  precursor  is  not  detected  easily  at  BING,  WND,  or  SANY, 
although  there  is  some  indication  of  anomalous  energy  preceding  DF  by  about  6-7  seconds  at 
DHN  (1 19.8°).  Figures  5  and  6  display  events  840806.1201  and  8501 1 1.1441,  respectively,  each 
illuminating  the  same  portion  of  the  P'  travel-time  curve  ( 1 32°  -  1 37°).  Precursive  arrivals  to  DF 
are  clear  in  Figure  5,  but  the  wavetrain  is  longer  (13  -  16  seconds)  and  more  emergent  than  in 
Figure  4,  with  amplitudes  increasing  toward  the  DF  arrival  (compare,  for  example,  WPNY  in 
Figures  4  and  5).  These  events  are  at  two  different  azimuths,  however,  and  it  is  not  yet  certain 
whether  we  are  observing  source  effects,  source  site  effects,  or  differences  in  the  characteristics  of 
D"  or  the  CMB  between  the  two  paths.  Observations  of  short-period  P  at  closer  distances  can 
assist  in  the  interpretation  of  this  difference.  We  note,  however,  that  the  onset  times  are 
qualitatively  consistent  with  the  theoretical  least  times  given  in  H addon  and  Cleary  \  1974,  Figure 
6].  We  contrast  840806.1201  in  Figure  5  with  850111.1441  in  Figure  6,  where  precursive 
arrivals,  if  present  at  all,  are  of  much  smaller  amplitude  relative  to  DF.  Again,  we  cannot  tel! 
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Figure  4:  Record  section  comprising  selected  recordings  on  NYSSA  of  event 
8601 15.2017,  h=  145  km.  DF-precursor  arrivals  arc  evident  as  distributed  waveforms  with 
impulsive  onsets  at  distances  beyond  122°.  Travel-time  curve  is  JB  with  constant  offset 
added  to  match  observed  DF  movcoul.  One  minute  of  record  is  shown. 


Figure  5:  Record  section  comprising  selected  recordings  ' 
840806.1201,  h=242  km.  Clear  prccursivc  arrivals  arc  evident  at 
the  onset  of  the  precursor  is  more  emergent  than  observed  in  Figur 
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Figure  6:  Record  section  comprising  selected  recordings  on  NYSSA  of  event 
8501 11  1441  h=189  km.  Though  having  nearly  the  same  range  and  back  azimuth  as  the 
record  section  shown  in  Figure  5,  the  DF-prccursors  shown  here  arc  much  less  evident. 
Note  also  that  the  P'  arrival  is  much  simpler. 
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Figure  7:  Record  section  comprising  selected  recordings  on  NYSSA  of  event 
850413.0106,  h=99  km.  The  large  amplitude  arrival  is  F'BC; DF  appears  as  the  much 
smaller  arrival  beyond  147°.  Travel-time  curves  arc  JB,  with  a  constant  offset  added  to 
match  the  DF  branch. 


without  further  study  whether  we  are  observing  the  effects  of  source  complexity  or  CMB 
heterogeneity  (the  DF  waveforms  in  Figure  6  are  much  simpler  than  those  in  Figure  5). 


I 

i 

!  Figure  7  displays  a  record  section  collected  just  beyond  the  fi-cusp  at  about  143°.  The 

dominant  arrival  is  P'BC,  with  appearing  only  as  a  small,  low-amplitude  arrival  beyond  145°. 
The  travel-time  curves  shown  are  calculated  for  the  JB  model,  with  a  constant  time  offset  of  2. 1  s 
added  in  order  to  line  up  P’DF  with  the  data.  The  P'D F  and  P'BC  travel  times  appear  to  be  well  fit  by 
the  offset  model  at  distances  greater  than  146°,  but  for  distances  closer  than  the  BC  -  DF  crossover, 
the  P'BC  arrivals  are  much  earlier  (about  2  seconds)  than  predicted. 

!  Figures  4  -  7  are  given  as  examples  of  the  data  quality  and  the  detail  available  when  station 

sites  are  densely  distributed  at  appropriate  epicentral  distances.  A  more  comprehensive  analysis  of 
these  data  is  underway  at  LDGO,  and  will  be  the  subject  of  another  publication. 

1  In  the  discussion  that  follows,  we  use  the  PKPDF  waveform  from  event  840806.1201  to 
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■  estimate  signal  coherences. 
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Figure  8:  Amplitude  coherence  between  two  different  sensors  situated  on  die  same  pier, 
for  one-minute  signal  duration,  using  the  PKP-DF  trigger  for  event  840806.1201. 
Cross-spcctrum  and  auto-spectra  estimates  arc  smoothed  with  a  .5  Hz  frequency  boxcar. 

Coherence  is  high  out  to  3  Hz,  and  acceptable  to  5  Hz.  Lower  values  of  coherence  at  low 
frequencies  arc  due  to  instrumental  drift. 

5.  Coherence  estimates 

For  a  period  of  time,  LDGO  operated  two  sensors  (Geospace  US- 10  and  a  Ranger,  both 
with  natural  periods  of  1  s)  on  the  same  pier  at  PAL,  using  amplifiers  and  filters  having 
characteristics  identical  to  the  rest  of  the  network.  Figure  8  shows  their  normalized  coherence 
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(smoothed  over  .5  Hz  spectral  window)  from  DC  to  10  Hz  for  the  PKPD F  signal  produced  by 
840806.1201.  The  coherence  is  near  unity  between  .5  and  2.5  Hz,  and  maintains  a  level  of  about 
.6  to  5  Hz.  There  is  some  incoherency  below  .5  Hz  which  is  noticeable  in  the  time  series  as 
long-period  drift.  This  is  easy  to  spot  upon  examination,  however,  and  is  not  a  factor  in  our 
interpretation  of  the  waveforms.  Evidently,  there  are  no  pathologic  instabilities  in  the  electronics 
which  might  impact  signal  coherence  below  2  Hz.  In  what  follows,  we  use  2  Hz  as  the  upper 
frequency  limit  for  our  analysis. 

We  concentrate  in  this  report  on  coherence  measurements  at  the  Adirondack  and  northern 
New  York  sites  (sub-nets  I  and  II  in  Figure  1),  using  six  vertical-component  stations  at  CTR, 
MSNY,  BGR,  WNY,  PNY,  and  SKT  (near  GNF).  For  the  most  part,  the  stations  are  situated  on 
hard-rock  with  some  glacial  till  cover  with  low  cultural  and  microseismic  noise  characteristics 
[LDGO  Field  reports].  Inter-station  spacing  ranges  from  24.7  km  to  128.8  km,  and  the  median 
spacing  is  97.7  km.  For  these  preliminary  studies,  we  extracted  one  minute  of  PKP  signal 
beginning  with  the  DF-precursor  and  tapered  the  time  series  with  a  10%  cosine.  We  assume  that 
the  instrument  responses  are  characterized  by  the  nominal  poles  and  zeros  given  in  Lee  and  Stewart 
[1981]. 
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Figure  9:  Inter-station  coherence  as  a  function  of  frequency  between  CTR  and  other 
stations  of  the  Adirondack  and  northern  New  York  subnets,  for  onc-minutc  signal 
duration.  Numbers  to  the  right  of  the  legend  give  inter-station  spacing  in  kilometers. 
Coherence  levels  arc  lower  for  increased  station  separation  as  expected. 


Figure  9  shows  coherence  estimates  (.5  Hz  spectral  window  smoothing)  between  CTR  and 
the  other  five  stations  as  a  function  of  frequency.  The  best  coherence  is  exhibited  for  CTR*SKT, 
which  also  have  the  closest  inter-station  spacing  (between  2  and  3  signal  wavelengths).  Even  so, 
the  coherence  drops  below  .5  above'  .75  Hz.  Not  surprisingly,  sub-array  heterogeneity  affects  the 


coherence  estimate.  Figure  10  displays  coherence  for  two  station  pairs  having  nearly  identical 
inter-station  spacing  (WNY*PNY  and  WNY*SKT).  Differences  in  the  coherence  curves  are 
indicators  of  differences  in  station  response  or  station  site  effects  (noise,  transmission  losses, 
differences  in  structure,  and  so  on).  The  curves  have  similar  fall-offs  above  about  .7  Hz,  though 
the  WNY*SKT  coherence  is  lower. 


840806.1201 


Figure  10:  Coherence  between  WNY,  PNY,  and  SKT.  WNY  is  54  km  from  both 
PNY  and  SKT.  PNY  is  109  km  from  SKT.  Differences  in  coherence  indicate  that  the 
distribution  of  heterogeneities  is  not  isotropic. 
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Figure  11:  Coherence  as  a  function  of  inter-station  spacing  for  each  of  several  narrow 
bands.  The  "hole"  between  50  and  60  km  is  probably  due  to  improper  instrument 
response. 


Figure  1 1  displays  coherence  as  a  function  of  inter-station  spacing  for  several  narrow  bands 


in  frequency.  The  coherence  "hole"  between  50  and  60  km  is  presumably  due  to  problems  with  the 
station  at  WNY  (1-second  vs.  2-second  seismometer,  site  characteristics),  though  nothing  in  the 
field  installation  notes  indicates  any  special  siting  problems  or  instrumental  defects.  Relatively  high 
coherence  is  seen  at  low-frequencies  (<.5  Hz)  out  to  distances  of  about  60  to  80  km  (3-5  crustal 
wavelengths),  somewhat  greater  than  observed  at  LASA  [Aki,  1973]  or  NORSAR  [Bungum  et  al., 
1971]. 

While  these  results  are  very  preliminary  (and  point  to  the  necessity  of  obtaining  appropriate 
instrument  response  functions),  they  nevertheless  show  the  potential  usefulness  of  a  well-situated 
array  for  teleseismic  structural  studies  and  subjacent  site  structure. 

6.  Plans  for  further  analysis 

Future  work  will  include  the  reconstruction  of  individual  station  calibration  responses 
(particularly  phase  responses)  from  field  installation  and  maintenance  notes.  We  will  also  develop 
more  robust  coherence  estimation  techniques,  and  include  more  events  in  our  coherence  estimates. 
The  one-minute  time  windows  used  in  our  estimation  include  not  only  primary  signal  but  some 
clearly  reverberatory  waveforms;  future  analyses  will  examine  shorter  time  windows  and  will  use 
polarization  analysis  to  detect  mode  conversions.  In  addition,  LDGO  is  installing  new  permanent 
instrumentation  and  acquiring  new  portable  instruments.  These  will  be  deployed  in  various 
configurations  for  both  intra-network  coherence  and  PKP  studies. 
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